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A  probabilistic  framework  for  prediction  of  material  property  distributions  from  small 
scale  (i.e.  2-grain)  models  is  proposed.  Monte  Carlo  Simulation  and  kernel  density 
estimation  are  used  to  estimate  the  material  property  distribution  of  a  grain  boundary  with 
a  2-grain  model.  Extreme  value  and  order  statistics  are  then  employed  to  estimate  the 
distribution  of  larger  microstructure  models.  An  example  of  the  methodology  is  presented 
for  identifying  the  applied  uniaxial  stress  at  which  plastic  slip  initiates  in  a  titanium  alloy 
with  a  crystal  elastic  finite  element  model.  The  framework  was  verified  by  comparing  the 
predicted  plastic  slip  initiation  strength  distribution  with  the  obtained  distribution  from 
Monte  Carlo  Simulation  of  larger  scale  finite  element  models  (i.e.  n-grain  models,  up  to  ~600 
grain  RVE).  The  methodology  performs  well  for  larger  microstructure  models  but  less  so  for 
smaller  ones,  for  a  much  smaller  computational  cost. 
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Nomenclature 

Applied  uniaxial  stress  necessary  for  plastic  slip  initiation 
slip  system  index 

Resolved  shear  stress  on  the  a th  slip  system 
Critical  resolved  shear  stress  of  the  CCth  slip  system 
Slip  plane  normal  unit  vector  of  the  ath  slip  system 
Slip  direction  unit  vector  of  the  Ctth  slip  system 

Cumulative  distribution  function  of  Op  of  the  2-grain  model  (parent  distribution) 

Number  of  grains  neighboring  a  grain  of  interest 

Cumulative  distribution  function  of  Op  of  the  weakest  pair  of  grains  in  a  collection  of  N  pairs 
Number  of  grains  in  a  model 

Cumulative  distribution  function  of  Op  of  the  weakest  grain  in  a  collection  of  n  grains 


I.  Introduction 

THE  focus  of  much  of  the  recent  research  in  computational  materials  has  been  in  digital  representation  of 
microstructures,  constitutive  relationships  of  different  microstructural  features,  and  multiscale  modeling.1  An 
important  aspect  of  multiscale  modeling  is  the  development  of  a  representative  volume  element  (RVE),  a  volume  of 
the  micro  structure  of  a  large  enough  size  that  the  distribution  of  the  response  of  interest  (i.e.  yield  strength)  is  in 
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agreement  with  the  distribution  of  the  macroscopic  material  response.  Determining  the  size  of  the  RVE  for  a  given 
response  has  been  the  focus  of  much  research.  Often  the  size  of  the  RVE  is  too  large  to  be  computationally  tractable 
for  property  identification.  For  this  reason,  statistical  volume  elements  (SVEs)  were  developed.  The  SVE  is  smaller 
in  general  than  the  RVE,  but  is  sampled  many  times  to  build  up  the  statistics  that  would  be  captured  by  the  RVE.2  A 
natural  question  that  arises  concerning  the  SVE  is:  How  small  can  an  SVE  be  with  minimal  effect  on  the  accuracy  of 
the  predicted  macroscopic  material  property  distribution?  The  smallest  an  SVE  (measured  in  number  of  grains) 
could  be  is  that  of  a  single  grain.  This  SVE  would  not  include  any  effects  on  the  response  due  to  grain  boundaries 
and  as  such  is  only  suitable  for  modeling  single  crystals.  The  next  simplest  microstructure  is  that  of  2  grains  -  a 
bicrystal.  For  some  macroscopic  material  properties  of  interest,  can  a  collection  of  2-grain  SVEs  be  used  to  build  up 
the  macroscopic  distributions?  This  work  describes  a  probabilistic  framework  that  attempts  to  obtain  such 
macroscopic  property  distributions  from  small  2  grain  SVEs. 

The  property  of  interest  in  this  study  is  the  plastic  initiation  strength  of  the  material.  For  a  given  applied  uniaxial 
load,  if  any  slip  systems  are  activated  at  any  location  in  the  material,  plastic  slip  has  initiated  at  that  location  and  the 
plastic  slip  initiation  strength  of  the  entire  model  has  been  reached.  This  type  of  behavior  is  called  Weakest  Link 
Theory  (WLT)  and  has  a  long  history  of  application  in  materials  science.  The  statistics  ofWLT  fall  into  the  category 
of  Extreme  Value  Theory  (EVT),  a  subset  of  Order  Statistics.  Initially  EVT  was  described  by  Weibull  for 
characterizing  the  rupture  of  solids3;  Batdorf  revisited  the  concepts  again  for  the  strength  of  ceramics4,  and  more 
recently  Tryon  and  Cruse  explored  these  concepts  for  fatigue  life5.  WLT  is  applied  often  in  systems  reliability 
analysis  for  series  systems.  The  approach  taken  here  in  this  work  is  to  view  the  microstructure  as  a  system  of  grains 
and  grain  boundaries,  each  of  which  have  some  probability  of  slip  initiation.  In  the  parlance  of  system  reliability,  for 
this  study  slip  initiation  corresponds  to  failure  of  the  material,  and  for  a  given  stress  state,  the  material  is  either  in  a 
state  of  failure  or  survival  (slip  or  no-slip).  A  concise  review  of  order  statistics  and  extreme  value  theory  is  given  is 
given  in  subsection  A.  The  model  for  slip  initiation  strength  of  a  titanium  alloy  is  presented  in  subsection  B. 

A.  Order  Statistics 

Order  statistics  is  concerned  with  the  distribution  of  the  rh  smallest  observation  in  a  sample  of  n  observations. 
For  a  given  sample  (X] ,  X2 , . . .  Xn ) ,  the  observations  may  be  sorted  in  increasing  order  and  relabeled 

(XVn ,  X2.n Xn.n  ).  The  rh  order  statistic  is  the  rth  smallest  value  and  labeled  Xrn  .  The  extreme  values,  the 
minimum  and  maximum,  are  then  XVn  and  X n.n  . 

If  the  sample  is  random,  then  uncertainty  in  X  is  described  by  its  distribution  Fp(x )  .  It  is  given  the  subscript  “P” 
to  identify  it  as  the  parent  distribution  that  determines  the  distributions  of  Xln , X2  . X  .  The  distribution  of 

the  rh  order  statistic  is  then  given  in  terms  of  Fp{x)  as 


=  2  V  *  "/F,w(r«n-r  +  1) 


(1) 


k=r' 


Extreme  Value  Theory,  a  subset  of  order  statistics,  gives  the  distribution  of  the  minimum  for  r  =  1  in  a  simpler  form 
as 


(2) 


B.  Model  of  Slip  Initiation  in  Ti-  a. 

Ti-6A1-4V  is  a  common  alloy  of  titanium  often  used  in  aircraft  structures  and  gas  turbine  engine  components.  Its 
high  strength,  low  density,  resistance  to  corrosion,  and  ability  to  withstand  higher  temperatures  has  lead  to  its 
widespread  application  '.  It  is  categorized  as  an  O.  +  /3  alloy,  and  contains  both  the  primary  alpha  phase  and  the  beta 
phase.  The  alpha  phase  has  a  hexagonal  close  packed  crystal  structure  (hep)  and  the  beta  phase  has  a  body  centered 
cubic  (bcc)  crystal  structure.  Depending  on  mechanical  processing  and  heat  treatment,  it  presents  in  a  microstructure 
ranging  from  fully  lamellar  to  bi-modal  to  fully  equiaxed,  with  varying  volume  fractions  of  a  and  /3.  The  specific 
microstructure  studied  here  is  fully  equiaxed.  The  beta  phase  was  not  modeled,  so  the  quantitative  results  of  this 
study  may  only  correspond  to  microstructures  with  very  high  volume  fractions  of  the  alpha  phase.  The 
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microstructure  modeled  is  a  collection  of  equiaxed  grains  made  of  Ti-  a.  Each  grain  has  a  different  crystallographic 
orientation,  as  seen  in  Figure  1. 


Figure  1.  Fully  Equiaxed  Ti-6A1-4V  Model  -  Each  grain  is 
assigned  a  different  crystallographic  orientation.  No  texture 
was  assumed. 


Each  grain  was  assigned  the  same  elastic  material  properties  in  the  form  of  a  stiffness  tensor.  The  hep  Ti-aP 
phase  is  modeled  as  transversely  isotropic.  The  elasticity  tensor  of  transversely  isotropic  materials  has  the  form: 
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where  Cn  =  C,2  and  C13  =  C23  and  C55  =  C66 ,  yielding  6  independent  constants  that  describe  the  elastic 
response  of  the  material.  The  elasticity  tensor  selected  for  this  study  was  that  given  in  [4]. 
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Table  1.  Anisotropic  Elastic  Constants  for  Ti-a  [GPa] 
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162.4 


180.7 


92.0 


69.0 


117 


35.2 


The  anisotropic  stiffness  and  varying  crystallographic  orientations  of  each  grain  give  rise  to  non-uniform  stresses 
in  the  microstructure,  even  in  the  presence  of  uniform  macroscopic  loading.  These  non-uniform  stresses  give  rise  to 
local  plastic  slip  in  the  material. 

The  hep  crystal  lattice  has  several  slip  system  families.  Slip  in  crystals  generally  occurs  on  planes  in  the  lattice 
with  relatively  high  close  packing  of  atoms.  These  planes  are  called  slip  planes  and  can  be  specified  by  the  vector 
normal  to  the  plane.  For  a  given  slip  plane,  the  slip  occurs  in  a  certain  slip  direction  when  the  shear  stress  resolved 
on  that  plane  in  that  direction  exceeds  the  critically  resolved  shear  stress  tCRSS  for  that  slip  system.  Slip  systems  that 
are  identical  under  a  rotation  of  the  crystal  lattice  due  to  crystallographic  symmetry  are  collected  in  a  set  called  a 
family.  The  three  slip  system  families  for  Ti-  ap  in  which  plastic  slip  is  most  commonly  observed  are  the  prismatic, 
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basal,  and  pyramidal.  These  slip  systems  for  Ti-  CCp  are  given  in  Miller-Bravais  indices  in  Table  2  and  presented  in 
Figure  3. 


Table  2.  Slip  Systems  for  Ti-aP 


Slip  System  Family 

Slip  Plane 

Slip  Direction 

Prismatic 

{lOlo}  (1120) 

(lloo) 

(ioTo) 

(olio) 

[ll20] 

[l210] 

[2110] 

Basal 

(0001) 

[1120] 

{0001}  (1120) 

(0001) 

[1210] 

(0001) 

[2110] 

(1101) 

[1120] 

(1 101) 

[1120] 

Pyramidal 

(1011) 

[12IO] 

{ion}  (1120) 

(Toil) 

[12IO] 

(01  Tl) 

[2110] 

(0111) 

[2110] 

Figure  2.  Prismatic,  Basal,  and  Pyramidal  Slip  System  Families  for  hep  Crystal  Lattices 

For  the  ath  slip  system,  the  resolved  shear  stress  on  the  slip  plane  with  unit  normal  n  and  slip  direction  b  is 
found  with 


<ss  =  h>ijnaj 


(4) 


The  slip  plane  normal  and  slip  direction  vectors  are  found  relative  to  the  local  material  coordinate  system.  They 
must  be  transformed  into  the  specimen  coordinate  system  to  be  implemented  with  the  stresses  from  ABAQUS  in  the 
above  equation. 

The  critically  resolved  shear  stress  for  each  slip  system  family  is  given  in  Table  3  below. 
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Table  3.  Critically  Resolved  Shear  Stress  for  Slip  System  Families  for  Ti-aP 


Slip  System  Family 

Tcrss  [MPa] 

Prismatic 

320 

Basal 

322 

Pyramidal 

846 

By  comparing  the  critically  resolved  shear  stress  with  the  resolved  shear  stress  given  an  applied  uniaxial  load, 
the  applied  uniaxial  load  necessary  for  plastic  slip  initiation  may  be  determined.  To  do  so,  a  performance  function 
for  each  slip  system  is  developed  such  that  when  its  value  is  less  than  or  equal  to  zero,  slip  has  initiated  on  that 
particular  system. 


8  ~^CRSS  T RSS  (5) 

g  <  0  :  Failure  (Plastic  Slip  Initiates) 
g  >  0  :  Survival  (No  Plastic  Slip) 


Assuming  a  linear  relationship  between  the  applied  load  and  the  value  of  the  resolved  shear  stress  for  each  slip 
system  at  every  location  in  the  micro  structure 


trss  =  ma+  b 


(6) 


The  unknowns  in  Eq.  7  are  m  and  b.  The  results  from  one  finite  element  analsyis,  TRSS  =  me 7  +  b  can  be  combined 
with  the  assumption  that  all  of  the  resolved  shear  stresses  are  zero  when  the  applied  uniaxial  stress  is  zero, 


0  =  m(0)  +  b.  It  is  then  evident  that  the  intercept  b  is  zero  and  the  slope  m  = 


The  extrapolation  for  the  slip 
O 


initiation  strength  as  shown  in  Figure  3  is  then 


_  T crss 

T RSS 


(7) 


This  Op  is  the  applied  uniaxial  load  necessary  to  initiate  plastic  slip.  A  value  for  Op  can  be  found  for  every 
location  in  the  model  and  the  location  with  the  smallest  value  of  Op  limits  the  allowable  uniaxial  load  that  can  be 
applied  to  the  model  without  initiating  slip.  This  weakest  link  behavior  suggests  application  of  extreme  value  theory. 
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Figure  3.  Extrapolation  for  the  Uniaxial  Stress 
Necessary  to  Initiate  Plastic  Slip  on  a  Given  Slip 
System  -  The  response  TRSS  is  obtained  from  a  single 
analysis  with  applied  load  o  ■  These  are  used  to 
extrapolate  for  Op  in  which  the  resolved  shear  stress 
equals  the  critical  resolved  shear  stress. 

The  computational  expense  of  RVE-size  models  prohibit  the  evaluation  of  large  sample  sizes  for  prediction  of 
material  property  distributions.  However,  a  2-grain  crystal  elastic  finite  element  model  is  small  enough  to  allow  the 
use  of  Monte  Carlo  Simulation.  The  uniaxial  stress  for  plastic  slip  initiation  may  be  efficiently  obtained  for  a  large 
number  of  realizations  to  estimate  the  parent  distribution.  The  extreme  value  theory  methodology  to  predict  the 
distribution  of  strength  of  larger  polycrystalline  specimens  is  described  in  Section  II.  The  results  from  the  titanium 
alloy  example  are  discussed  in  Section  III  and  the  conclusions  are  presented  in  Section  IV.  Despite  the  narrow  scope 
involved  with  modeling  only  the  alpha  grains  in  a  fully  equiaxed  microstructure,  the  extreme  value  theory 
framework  for  identifying  material  properties  will  be  demonstrated.  The  framework  can  be  further  extended  to  other 
microstructures  with  appropriate  modifications. 

II.  Methodology 

The  probabilistic  framework  developed  in  this  work  builds  the  distribution  of  the  response  of  the  macroscopic 
material  from  the  distribution  of  the  response  of  an  SVE  using  order  statistics.  The  key  steps  in  the  methodology  are: 

A.  Sample  the  small-scale  model  with  Monte  Carlo  simulation  (MCS)  to  develop  the  distribution  of 
the  SVE  response 

B.  Apply  Extreme  Value  Theory  and  Order  Statistics  to  estimate  the  distribution  of  the  response  of 
the  large-scale  model,  the  RVE,  and  hence  the  macroscopic  material,  from  the  SVE  response 
distribution 

Each  of  these  steps  are  discussed  in  detail  in  the  following  subsections. 

A.  MCS  of  small-scale  model  for  SVE  response  distribution 

A  small-scale  finite  element  model  was  developed  with  the  ABAQUS  software  and  in  this  work  is  called  an 
SVE.  In  this  preliminary  study,  the  SVE  is  a  2  grain  or  bicrystal  microstructure.  Each  grain  is  modeled  as  a  cube, 
and  their  shared  face  makes  up  their  grain  boundary.  A  perfect  bond  was  assumed  at  the  grain  boundary.  The  pair  of 
grains  were  enveloped  in  a  third  grain  that  was  given  the  isotropic  macroscopic  material  properties  of  alpha  titanium 
to  approximate  the  constraints  on  the  grains  that  would  occur  in  a  full  micro  structure.  The  model  was  meshed  with  a 
grid  of  20x30x40  elements.  The  geometry  is  presented  below  in  Figure  4. 
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Figure  4.  SVE  Geometry  -  The  anisotropic  grains  are  red  and  blue, 
the  isotropic  material  is  gray. 


Monte  Carlo  simulation  (MCS)  is  performed  with  the  small-scale  model  and  the  distributions  of  the  input 
variables:  the  crystallographic  orientation  of  each  grain.  For  each  sample,  a  random  rotation  was  applied  to  the 
material  orientation  of  each  grain  to  generate  that  grain’s  crystallographic  orientation  using  unit  quaternions6.  The 
sampling  routine  is  included  here  for  completeness. 

1.  Sample  three  uniform  random  variables:  Uv  U2,  W3~U[0,1] 

2.  Compute  the  quaternion:  q  =  (/?0  ,Pi  ,/?2  ,/?3  ) 

P0  =  -  uA  ■  sin(2 jvu2) 


f\  =  -Jl  -  M|  ■  cos(2 jiu2) 


P2  =  sin(2jTM3) 


/?3  =  cos(2 jtu3) 


(8a,b,c,d) 


3.  Compute  the  associated  rotation  matrix,  R  that  transforms  a  coordinate  in  the  local  material  coordinate 
system  to  one  in  the  global  specimen  coordinate  system. 


R  = 


Pi +  -  Pi  -  Pi  +  2(m-f\A) 

2 (A A  -  A.A )  -  A2  +  A2  -  A2  2<7bA  +  PA) 

2(A  A  +  PA)  2 (P2p3  -  PA )  Pi  -  Pi  -  Pi  +  Pi 


(9) 


4.  The  unit  vectors  that  define  the  coordinate  system  of  the  material  in  the  global  coordinate  system  are  then: 
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Xsx 

1 
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Xsy 
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0 
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ysx 

0 
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-=m 

1 
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0 

7 

^sx 

0 

{usz}  =  [R]{umz}^< 

7 

^ sy 

=  [Rl 

0  > 

7 

l  SZ  J 

1 

These  vectors  are  used  in  ABAQUS  to  assign  the  material  coordinate  systems.  The  next  step  was  to  obtain  the 
response  of  the  model  for  each  sample. 

The  linear  elastic  finite  element  analysis  was  performed  for  each  sample  to  obtain  the  stress  state,  O^-,  at  each 
integration  point  in  the  model.  It  is  assumed  that  plastic  slip  would  first  initiate  at  the  grain  boundary  because  of  the 
discontinuity  in  the  material  properties,  and  thus  only  the  integration  points  in  elements  along  the  grain  boundary 
were  considered.  The  stess  tensor  of  the  integration  points  along  the  grain  boundary  were  used  to  calculate  the 
resolved  shear  stress  on  each  of  the  slip  systems  for  primary  alpha  titanium.  This  was  then  compared  with  the 
critically  resolved  shear  stress  for  each  slip  system.  An  example  of  the  stresses  is  given  in  Figure  5. 


Figure  5.  Example  SVE  Stresses  gxx  -  Only  the  stresses  along  the 
grain  boundary  are  used  in  the  slip  initiation  analysis.  A  uniaxial 
stress  was  applied  on  and  normal  to  the  near  face  on  the  right. 
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The  slip  initiation  strength  was  taken  as  the  smallest  uniaxial  load  necessary  to  activate  any  slip  system  at  any 
integration  point  in  the  model.  With  a  large  sample  size,  kernel  density  estimation  is  then  used  to  obtain  the 
distribution  of  the  slip  initiation  strength  of  a  bicrystal. 

B.  Estimate  RVE  response  distribution  with  Order  Statistics  and  SVE  response  distribution 

Development  of  the  RVE  response  distribution  is  facilitated  by  the  application  of  order  statistics  and  the  SVE 
response  distribution. 

The  distribution  of  the  slip  initiation  strength  of  a  bicrystal  is  obtained  in  the  previous  step.  This  2-grain 
distribution  is  then  used  as  an  approximation  of  the  slip  initiation  strength  for  an  arbitrary  pair  of  neighboring  grains 
in  the  bulk  material.  The  key  assumption  made  in  this  methodology  is  that  the  bicrystal  distribution  is  a  good 
approximation  of  the  distribution  of  the  behavior  of  neighboring  grains  in  the  bulk  material. 

Based  on  this  approximation,  the  SVE  distribution  is  then  the  parent  distribution  for  use  with  extreme  value  and 
order  statistics  to  develop  the  distribution  of  the  slip  initiation  strength  of  the  RVE.  The  parent  distribution  is 
denoted  by  its  CDF,  Fp  (o  p  ) . 

The  distributions  for  order  statistics  and  EVT  given  in  equations  1  and  2  are  for  the  case  where  the  n  samples  are 
independently  and  identically  distributed.  Applying  these  equations  to  the  distribution  of  the  initiation  strength  of 
different  grains  in  an  SVE  will  introduce  some  bias  because  the  strength  distribution  of  one  grain  depends  on  the 
same  factors  that  also  effect  other  neighboring  grains  (i.e.  Schmid  Factor)  and  their  distributions  are  not 
independent.  It  is  not  clear  how  much  bias  is  introduced  by  making  the  independence  approximation.  The 
verification  tests  may  shed  light  on  this  issue,  and  even  indicate  if  the  proposed  methodology  can  at  least  provide  a 
bound  on  the  response  distribution. 

1.  Application  of  Order  Statistics  to  determine  the  RVE  slip  initiation  strength  distribution  from  the  SVE  slip 
initiation  strength  distribution 

Given  the  parent  distribution  (i.e.  the  SVE  slip  initiation  strength  distribution  CDF),  Fp(<Jp )  and  extreme  value 
theory,  the  distribution  of  the  slip  initiation  strength  of  a  single  grain  with  N  adjacent  grain  pairs  as  the  minimum 
strength  of  the  N  pairs  and  modifying  Eqn.  2. 

FiM(oP)-l-(l-FP(oPjf  (ll) 

Eqn.  11  assumes  the  value  of  N  is  deterministic.  Given  some  grain  in  the  material,  the  number  of  its  neighboring 
grains,  N  is  a  random  variable  with  a  probability  mass  function*.  An  example  of  one  is  given  in  Figure  ##. 


Figure  6.  PMF  of  Number  of  Nearest  Neighbors  -  Uncertainty  in  the 
number  of  neighboring  grains  for  a  given  grain  of  interest  can  be 
accounted  for  with  the  total  probability  theorem. 
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The  uncertainty  in  N  may  be  addressed  with  the  total  probability  theorem*.  The  total  probability  theorem  is  given  by 

p(a)-p(a\e,)p(e, )  +  ...  + p(a\en)p(e„)-'2p(a\e,)p(e,)  (12) 

N 

Eqn.  11  is  modified  with  Eqn.  12  to  account  for  the  uncertainty  in  the  number  of  neighboring  grains  of  a  randomly 
selected  grain  with 


max  N  -  i 

Fl:Jv(°p)=  2  (13) 

N= 1  L 

where  PN  is  the  probability  that  a  grain  will  have  N  neighbors.  This  gives  a  model  of  the  distribution  of  the  slip 
initiation  strength  of  a  single  grain. 

Using  order  statistics,  the  applied  uniaxial  stress  such  that  r  independent  grains  in  an  n  grain  model  experience 
slip  initiation  is 


(14) 


where  /  is  the  incomplete  beta  function. 


III.  Results  and  Discussion 

The  distributions  for  the  slip  initiation  strength  predicted  with  the  EVT  framework  for  various  size  SVEs  are 
presented  and  compared  with  results  from  MCS  of  SVE  FE  analyses  . 

A.  Distributions  predicted  by  2-grain  model 

The  2-grain  model  was  exercised  with  Monte  Carlo  Simulation  for  200  samples.  A  kernel  density  estimate  based 
on  this  data  was  used  to  obtain  Fp(op )  .  The  corresponding  PDF  and  CDF  are  displayed  in  Figure  7  below. 
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Slip  Initiation  Strength  [MPa] 


Figure  7.  PDF  and  CDF  of  Slip  Initiation  Strength  for  a  Random 
Pair  of  Grains  Predicted  by  2-grain  Model  -  It  is  assumed  the  2- 
grain  model  is  a  good  approximation  of  a  pair  of  grains  in  a 
polycrystal.  This  “parent”  distribution  is  used  to  predict  the 
distribution  for  a  polycrystal. 

The  2-grain  parent  distribution  was  was  then  used  in  Eq.  1 3  with  the  distribution  for  the  nearest  neighbors  given 
in  Figure  6  to  obtain  the  distribution  for  the  strength  associated  with  the  boundaries  about  a  random  grain.  The  PDF 
and  CDF  are  given  in  Figure  8. 


Figure  8.  PDF  and  CDF  of  Slip  Initiation  Strength  for  a  Random 
Grain  Predicted  by  2-grain  Model  -  This  is  the  “child”  distribution 
predicted  by  the  2-grain  model.  It  is  used  as  a  “ parent  ”  distribution  to 
predict  the  distribution  for  a  polycrystal. 

The  predicted  grain  strength  distribution  was  then  used  as  a  parent  distribution  in  Eq.  14  to  obtain  the 
distribution  of  the  slip  initiation  strength  of  a  polycrystal,  seen  in  Figure  9. 
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RVE  Strength 


Figure  9.  PDF  and  CDF  of  Slip  Initiation  Strength  for  a  Polycrystal 
Predicted  by  2-grain  Model  -  This  is  the  “child”  distribution  of  the 
child  predicted  by  the  2-grain  model.  It  was  developed  for  a  polycrystal 
of  30  grains. 


Slip  Initiation  Strength  [MPa] 


Figure  10.  Comparison  of  Predicted  Slip  Initiation  Strength 
Distributions  -  The  mean  and  variance  decrease  for  each  successive 
level.  The  predicted  RVE  distribution  is  for  a  polycrystal  of  30-grains. 


All  of  the  2-grain  model-predicted  slip  initiation  strength  distributions  have  a  mean  less  than  the  macroscopic 
yield  strength  for  titanium,  as  expected.  A  stronger  verification  of  the  framework  was  performed  by  running  crystal 
elastic  finite  element  analyses  on  many-grained  models  for  comparison  as  described  in  the  next  section. 


B.  Distributions  of  larger  SVEs 

Verification  tests  with  realizations  of  SVEs  with  more  than  2  grains  were  performed  with  ABAQUS. 
Realizations  of  the  statistically  equivalent  microstructures  were  generated  in  voxel  form  with  DREAM. 3D.  The 
voxels  were  then  used  as  brick  elements  in  a  finite  element  mesh  for  analysis  with  ABAQUS.  Three  different 
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volumes  were  considered  to  assess  the  2-grain  model’s  prediction:  0.01x0.01x0.01  mm,  0.03x0.03x0.03  mm,  and 
0.1x0. 1x0.1  mm. 

1.  ~6-grain  SVE 

The  first  verification  example  is  shown  in  Figure  11a.  It  is  a  cube  with  side  of  length  0.01mm.  The  number  of 
realizations  generated  for  the  verification  test  was  1,250.  The  number  of  grains  in  each  realization  varied  as  seen  in 
Figure  11c.  The  mean  number  of  grains  was  5.6  grains  and  the  standard  deviation  was  1.8  grains.  A  correlation  of 
-0.28  was  found  between  the  extrapolated  slip  initiation  strength  and  the  number  of  grains  in  the  SVE  as  seen  in 
the  scatter  plot  of  1  Id. 


Figure  lla,b,c,&d.  Example  0.01x0.01x0. 01mm  SVE  -  A  uniaxial  tensile 
stress  was  applied  on  and  normal  to  the  front  right  face  and  the  opposite  face 
(not  visible)  was  held  fixed.  The  slip  initiation  site  is  enclosed  by  the  red  circle 
and  occurred  at  a  grain  boundary’  for  this  geometry.  The  number  of  grains  in 
each  realization  was  random,  and  a  small  correlation  (-0.28)  was  found 
between  slip  initiation  strength  and  number  of  grains  in  the  SVE. 


The  obtained  distribution  of  the  slip  initiation  strength  for  the  ~6-grain  model  from  MCS  and  is  presented  in 
Figure  12a.  The  mean  values  was  900  MPa  and  the  standard  deviation  was  901  MPa.  The  MCS  distribution  is 
compared  to  the  predicted  distribution  from  EVT  in  Figure  12b.  The  predicted  distribution  has  a  mean  of  435  MPa 
and  standard  deviation  of  94  MPa.  The  predicted  distribution  differs  significantly  from  the  actual.  The  lower  bounds 
of  each  is  similar  and  the  predicted  distribution  is  conservative,  but  the  means  differ  by  465  MPa  and  the  EVT 
predicted  variance  is  not  conservative. 
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Figure  12a.  MCS  Distribution  of  Plastic  Slip  Initiation  Strength  of 
0.01x0. 01x0.01mm  SVE  -  Sample  size  =  1250  realizations,  mean  =  900  MPa, 
standard  deviation  =  901  MPa.  b.  Comparison  with  EVT  Predicted 
Distribution  -  mean  =  435  MPa,  standard  deviation  =  94  MPa.  The  two 
distributions  differ  significantly. 


2.  ~38-grain  SVE 

The  second  verification  example  is  presented  in  Figure  13a  and  is  a  cube  with  side  of  length  0.03mm.  The 
number  of  realizations  generated  for  the  verification  test  was  1,185.  The  number  of  grains  in  each  realization  varied 
as  seen  in  Figure  13c.  The  mean  number  of  grains  was  38  grains  and  the  standard  deviation  was  5.1  grains.  A 
negligible  correlation  of  —0.085  was  found  between  the  extrapolated  slip  initiation  strength  and  the  number  of 
grains  in  the  SVE  as  seen  in  the  scatter  plot  of  13d. 


Figure  13a,b,c,&d.  Example  0.03x0.03x0. 03mm  SVE  -  A  uniaxial  tensile 
stress  was  applied  on  and  normal  to  the  front  right  face  and  the  opposite  face 
(not  visible)  was  held  fixed.  The  slip  initiation  site  is  enclosed  by  the  red  circle 
and  occurred  at  a  grain  boundary’  for  this  geometry.  The  number  of  grains  in 
each  realization  was  random,  and  a  neglible  correlation  (-0.085)  was  found 
between  slip  initiation  strength  and  number  of  grains  in  the  SVE. 
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The  obtained  distribution  of  the  slip  initiation  strength  for  the  ~38-grain  model  from  MCS  and  is  presented  in 
Figure  14a.  The  mean  value  was  442  MPa  and  the  standard  deviation  was  61.5  MPa.  The  MCS  distribution  is 
compared  to  the  predicted  distribution  from  EVT  in  Figure  14b.  The  predicted  distribution  has  a  mean  of  360  MPa 
and  standard  deviation  of  37  MPa.  The  predicted  distribution  differs  slightly  from  the  actual.  The  lower  bounds  of 
each  is  similar  and  the  predicted  distribution  is  conservative,  and  the  means  differ  by  80  MPa.  However,  the  EVT 
predicted  variance  is  not  conservative. 


Figure  14a.  MCS  Distribution  of  Plastic  Slip  Initiation  Strength  of 
0.03x0. 03x0. 03mm  SVE  -  Sample  size  =  1185  realizations,  mean  =  442  MPa, 
standard  deviation  =  62  MPa.  b.  Comparison  with  EVT  Predicted 
Distribution  -  mean  =  360  MPa,  standard  deviation  =  37  MPa.  The  two 
distributions  differ  slightly. 


3.  ~62 5-grain  RVE 

At  the  time  of  this  writing,  only  1  full  scale  model  was  fully  analyzed  for  slip  initiation  strength.  The  equiaxed 
microstructure  is  shown  in  Figure  15.  It  is  a  cube  with  side  of  length  0.1mm  and  contains  625  grains.  The  uniaxial 
stress  necessary  to  initiate  plastic  slip  was  found  to  be  320.5  MPa,  with  slip  occurring  on  the  basal  plane.  The  EVT 
predicted  distribution  has  a  mean  of  290  MPa  and  a  standard  deviation  of  73  MPa.  The  mean  differs  from  the  sample 
point  by  30.5  MPa,  well  within  1  predicted  standard  deviation. 


Figure  15a.  Largest  Model  Examined,  0.1x0.1x0.1mm  RVE  -  A  uniaxial  tensile  stress  was  applied  on 
and  normal  to  the  front  right  face  and  the  opposite  face  (not  visible)  was  held  fixed.  Sample  size  =  1 
realization,  Slip  Initiation  Strength  =  320.5  MPa.  b.  Comparison  with  EVT  Predicted  Slip  Initiation 
Strength  Distribution  -  mean  =  290  MPa,  standard  deviation  =  73  MPa.  The  EVT  predicted  distribution 
encompasses  the  sample  point  of  the  large  scale  model. 
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As  discussed  previously,  it  is  expected  that  local,  microscopic  yielding  will  take  place  at  a  lower  applied  uniaxial 
stress  than  that  of  macroscopic  yielding.  The  range  of  the  macroscopic  yield  strength  of  Ti-6A1-4V  is  -800-1000 
MPa,  so  the  results  for  the  large  scale  model  fall  well  below  the  upper  limit  imposed  by  the  macroscopic  yield 
strength. 


IV.  Conclusion 

In  summary,  an  extreme  value  theory  framework  was  presented  for  predicting  polycrystalline  model  material 
property  distributions  from  bicrystal  model  experiments.  The  methodology  performs  poorly  for  predicting  the 
distribution  of  slip  initiation  strength  of  relatively  small  scale  micro  structure  models,  yet  performs  well  for  larger 
microstructure  models.  Potential  causes  of  discrepancy  between  the  EVT  predicted  distribution  and  the  actual 
distribution  may  be  due  to  non-local  effects  or  the  fact  that  grain  size  and  shape  are  ignored  as  of  this  writing  in  the 
2-grain  model.  Incorporating  their  effects  may  improve  the  ability  of  the  proposed  framework  to  predict  material 
property  distributions. 

The  successful  demonstration  of  the  proposed  framework  for  the  larger  models  indicates  the  methodology  may 
serve  as  a  useful  technique  for  estimating  certain  material  properties  at  a  much  smaller  computational  cost  than  that 
of  exercising  full  large  scale  models. 
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